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STABILITY OF A FINITE VOLUME SCHEME 
FOR THE INCOMPRESSIBLE FLUIDS 

Sebastien Zimmermann 1 

Abstract. We introduce a finite volume scheme for the two-dimensional incompressible Navier-Stokes 
equations. We use a triangular mesh. The unknowns for the velocity and pressure are respectively 
piecewise constant and affine. We use a projection method to deal with the incompressibility constraint. 
We show that the differential operators in the Navier-Stokes equations and their discrete counterparts 
share similar properties. In particular we state an inf-sup (Babuska-Brezzi) condition. Using these 
properties we infer the stability of the scheme. 

Resume. Nous introduisons ici un schema volumes finis pour les equations de Navier-Stokes in- 
compressibles en deux dimensions. Les maillages consideres sont formes de triangles. Les inconnues 
associees a la vitesse et la pression sont respectivement constantes et affines par morceaux. Nous 
utilisons une methode de projection pour traiter la contrainte d'incompressibilite. Nous verifions que 
les operateurs differentiels apparaissant dans les equations de Navier-Stokes et leurs analogues dis- 
crets verifient des proprietes similaires. Nous prouvons en particulier une condition inf-sup. Nous en 
deduisons la stabilite du schema. 
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1. Introduction 

We consider the flow of an incompressible fluid in a polyhedral set ft C R 2 during the time interval [0, T\. 
The velocity field u : f2 x [0, T] — > R 2 and the pressure field p : f2 x [0, T] — > R satisfy the Navier-Stokes equations 

u t --J-Au + (u-V)u + Vp = f, (1.1) 
Ke 

divu = 0, (1.2) 

with the boundary and initial conditions 

u| 9 n=0, u| t=0 =u . 

The terms Au and (u- V)u are associated with the physical phenomena of diffusion and convection, respectively. 
The Reynolds number Re measures the influence of convection in the flow. For equations (|1.1[) - (|1.2[) , finite 
clement and finite difference methods are well known and mathematical studies are available (see [9] for example) . 



Keywords and phrases: Incompressible fluids, Navier-Stokes equations, projection methods, finite volume. 
17 rue Barreme - 69006 LYON, e-mail: Sebastien.Zimmermann@ec-lyon.fr 

© EDP Sciences, SMAI 1999 



2 



TITLE WILL BE SET BY THE PUBLISHER 



For finite volume schemes, numerous computations have been conducted ( [12] and [1] for example). However, 
few mathematical results are available in this case. Let us cite Eymard and Herbin [6] and Eymard, 
Latche AND Herbin [7]. In order to deal with the incompressibility constraint (|1.2[) , these works use a 
penalization method. Another way is to use the projection methods which have been introduced by Chorin [4] 
and Temam [13]. This is the case in Faure [8] where the mesh is made of squares. In Zimmermann [14] the 
mesh is made of triangles, so that more complex geometries can be considered. In the present paper the mesh 
is also made of triangles, but we consider a different discretization for the pressure. It leads to a linear system 
with a better-conditioned matrix. The layout of the article is the following. We first introduce in section [2] 
the discrete setting. We state (section I2.1j) some notations and hypotheses on the mesh. We define (section 
12. 2| the spaces we use to approximate the velocity and pressure. We define also (section I2.3j) the operators 
we use to approximate the differential operators in (|1.1[) - (|1.2[) . Combining this with a projection method, we 
build the scheme in section [3] In order to provide a mathematical analysis, we show in section [4] that the 
differential operators in (|l.l[) - (|1.2p and their discrete counterparts share similar properties. In particular, the 
discrete operators for the gradient and the divergence are adjoint. The discrete operator for the convection term 
is positive, stable and consistent. The discrete operator for the divergence satisfy an inf-sup (Babuska-Brezzi) 
condition. From these properties we deduce in section [5] the stability of the scheme. 

We conclude with some notations. The spaces (L 2 , |.|) and (L°°, \\-\\oo) are the usual Lebesgue spaces and we 
set Lq = {q £ L 2 ; J n q(x)dx = 0}. Their vectorial counterparts are (L 2 , |.|) and (L°°, ||.[[oo) with L 2 = (L 2 ) 2 
and L°° = (L°°). For k £ N*, {H k , || • \\/.) is the usual Sobolev space. Its vectorial counterpart is (H fc , ||.||fe) with 
H fe = (H k ) 2 . For k = 1, the functions of H 1 with a null trace on the boundary form the space Hj. Also, we set 
Vu = (Vmi, Vu 2 ) t if u = (ui, u 2 ) G H 1 . If X C L 2 is a Banach space, we define C(0, T; X) (resp. L 2 (0, T; X)) 
as the set of the applications g : [0,T] — > X such that t — > |g(t)| is continuous (resp. square integrable). The 
norm ||.||c(o,T;X) is defined by ||g||c(o,T;X) = su Psg[o,t] ls( s )l- I n au calculations, C is a generic positive constant, 
depending only on fl, uq and f . 



2. Discrete setting 

First, we introduce the spaces and the operators needed to build the scheme. 

2.1. The mesh 

Let Th be a triangular mesh of Q. The circumscribed circle of a triangle K £ Th is centered at Xi<- and has 
the diameter hx- We set h = msx.KoX h hx- We assume that all the interior angles of the triangles of the mesh 
are less than ~, so that x#- 6 K. The set of the edges of the triangle K £ Ti is 8k ■ The symbol n.K,a denotes 
the unit vector normal to an edge a £ £k and pointing outward K. We denote by £h the set of the edges of the 
mesh. We distinguish the subset £ l h nt C £h (resp. Eff 1 ) of the edges located inside f2 (resp. on dQ). The middle 
of an edge a £ Eh is x CT and its length |c|. For each edge a £ £ h nt , let K a and L a be the two triangles having 
a in common. We set d a — d(x^,XL„). For all a £ £h Xt , only the triangle K a located inside is defined and 
we set d a = c^x^jXg.). Then for all a £ Eh we set r a = -p-. As in [5] we assume the following on the mesh: 
there exists C > such that 

VaeSh, d a >C\a\ and \a\ > C h. 

It implies that there exists C > such that 

Vcreff*, r a = \a\jd a > C. (2.1) 

2.2. The discrete spaces 

We first define 

P = {q £ L 2 ; Vif G %, q\x is a constant} , P = (Po) 2 - 
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For the sake of concision, we set for all qh G Po (rcsp. Vh G Po) and all triangle K G 7^: qk = qh\n (resp. 
vk = v/j|k)- Although Po ^ H 1 , we define the discrete equivalent of a H 1 norm as follows. For all Vh G Po 
we set ^ ^ 

||v fc || h = ( Yl t„\v l „-vkJ 2 + Yl t °M 2 I • (2-2) 
We have [5] a Poincare-like inequality: there exists C > such that for all v/j G Po 

|v fc |<C||v fc ||fc. (2.3) 
We also have [14] an inverse inequality: there exists C > such that for all Vh G Po 

ft||vfc|| h <C|vfc|. (2.4) 
From the norm ||.||^ we deduce a dual norm. For all v h G Po we set 

KIU,^ sup ^ftA. (2.5) 

For all u/j G Po and G Po we have (uh,Vh) < ||ufc||— ||v/j||k. We define the projection operator IIp : 
L 2 — > Po as follows. For all w G L 2 , ITp w G Po is given by 

VKeT h , (n Po w)| K = ^^w(x)dx. (2.6) 

We easily check that for all w G L 2 and v/j G Po we have (IIp w, -v h ) = (w, v^). We deduce from this that IIp 
is stable for the L 2 norm. We define also the operator n Po : H 2 — ► P . For all w G H 2 , IIp w G Po is given by 

VK€T h , ff Po w| K = w(x x ). 

According to the Sobolev embedding theorem, w G H 2 is a.e. equal to a continuous function. Therefore the 
definition above makes sense. We introduce also the finite element spaces 

P* = {v G L 2 ; ViC G %, v\ K is affine} , 
Pr = {v h G P* ; Va G 4"*, v h \ K „{*„) = v h \ L M , 
Pi = {vft G (Pi) 2 ; Vh is continuous and Vh\an = 0}. 

We have PJ C Hj. We define n P; : Hj -> PJ. For all v = («i, v 2 ) G Hj, IIpcv = w 2 ) G Pf is given by 

2 2 

= (^,$) G PJ , ]T (Vt£, V^) = E (V«i, V^). 

i=l i=l 

The operator Upc is stable for the H 1 norm. One checks ( [2] p. 110) that there exists C > such that for all 
v G H 1 

|v-n P cv| < Cft||v||i. (2.7) 

Let us address now the space P" c . If qh G P" c , we have usually Vqu G - L 2 . Thus we define the operator 
Vh ■ Pi c -^> Po by setting for all qh G Po and all triangle K &Th 

V h qh\ K = j^ jf_V(fodx. (2 ' 8) 



4 



TITLE WILL BE SET BY THE PUBLISHER 



The associated norm is denned by 

\Wh\\x,h = {\<lh\ 2 + \V h q h \ 2 ) 1 ' 2 . 
We have a Poincare-like inequality : there exists C > such that for all qu E P" c H Lq 

\qh\<C\V h q h \. (2.9) 

We define the projection operator Hpnc, For all q E H , Hpncq is given by 

Ver G £h , / (Upneq) da = qda. 

J (j J a 

One checks ( [2] p. 110) that there exists C > such that 

\ P -iip r p\<ch\\p\\ u IvfcCp-n^-p)! < (2.io) 

Finally, we use the Raviart-Thomas spaces (see [3]) 

RTq = {v h 6 Pi 1 ; \/<j<eSk, Vfc,|# • xi k ,<t is a constant, and v ft • n| 9f2 = 0} , 
RT = {v ft eRTg; VK G T h , V<r G £ x . v h | Ktr • n^ !<T = v A |i, • njf, ><7 }. 

For all G RT , A' G ?h and cr G £k we set (v^ • n# !<7 ) ff = v^|^ • n# i(7 . We define the operator Urt : H 1 — > 
RT . For all v G H 1 , n RTo v G RT is given by 



MK&T h , VcteEk, (n RTo v • n K , a )a = i— f / vdcr - 



(2.11) 



2.3. The discrete operators 

The equations (|l.l[l - (|1.2[l use the differential operators gradient, divergence and laplacian. Using the spaces 
of section |2"T21 we define their discrete counterparts. The discrete gradient : P™ c — > Po is defined by (|2.8p . 
The discrete divergence operator div/j : P — > P™ c is built so that it is adjoint to the operator V^. We set for 
all Vh G Po and all triangle K G Th 

3 I fX I 

G 4"*= ( div ft v /l )(x (T ) = — ^— (v Lct - v K J • njc )<F ; 

VaeEf?, (div h v h )(x a ) = - ]K V°\ L , vk„ ■ n V . (2.12) 

The first discrete laplacian : P" c — > P™ c ensures that the incompressibility constraint (|1.2p is satisfied in a 
discrete sense (see the proof of proposition 13.11 below) . We set for all q^ G P" c 

A h q h = div h (V h q h ). 

The second discrete laplacian A^ : Po — > Po is the usual operator in finite volume schemes [5]. We set for all 
V/i E Pq and all triangle K £ Th 
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In order to approximate the term (u • V)u in we define a bilinear form b/j : RTo x Po — » Po using the 
well-known upwind scheme [5]. For all % £ P , vj, £ P , and all triangle K £ Th we set 

b h (u h> v h )\ K = ^— l°1 (( u - n K,<r)t Vk+ (u-n Kj(T )- v L „y (2.13) 

We have set a + = max(a, 0), a~ = min(a, 0) for all a £ K. Lastly, we define the trilinear form b^ : RT x P x 
Po — > R 2 as follows. For all u/j 6 RT , V/j G Po, £ Po, we set 

b/^u^v/^w/j) = Y \K\w K - b h (u h ,v h )\ K . (2.14) 



3. The scheme 



In order to deal with the incompressibility constraint (|1.2p we use a projection method. This kind of method 
has been introduced by Chorin [4] and Temam [13]. The basic idea is the following. The time interval [0, T] 
is split with a time step k: [0,T] = \Jn =0 [t n ,t n+ i] with N £ W and t n = nk for all n £ {0, . . . ,N}. For all 
m £ {2, . . . , N}, we compute (see equation (|3.2p below) a first velocity field u™ ~ u(t TO ) using only equation 
(jl.lj) . We use a second-order BDF scheme for the discretization in time. We then project u™ (see equation (|3.4[) 
below) over a subspace of Po- We get a a pressure field ~ p{t m ) and a second velocity field u™ ~ u(i m ), 
which fulfills the incompressibility constraint (|1.2|) in a discrete sense. The algorithm goes as follows. For all 
m £ {0, . . . , N}, we set f™ = IIp f (t m ). Since the operator IIp is stable for the L 2 -norm we get 



|fn - |n Po f(t m )| < |f(t m )l < llf|l c(0 , T; L 2 )- ( 3J ) 

We start with the initial values 



u£eP nRT , uJeP„nRT pJeF n4 

For all n £ {1, . . . , N}, (u£ +1 ,p£ +1 , u£ +1 ) is deduced from (u£,p£, u%) as follows. 

• uJJ +1 <E Po is given by 

o yi+l _ 4 ri , n-l -, 

2^ A ^ u ft + b ^ u (,- u k i u h ) + VfcPft - t h , (3.2) 

• £ P? c n Lq is the solution of 

A.(K +1 ~K) = ^div,< +1 , (3.3) 



• u^ +1 e P is deduced by 



u „+l = -„+l _^ Vft(p n+l _ p ny (3 4) 

Existence and unicity of a solution to equation ()3.2j) is classical ( [5] for example) . The convection term in (|3.2p 
is well defined thanks to the following result. 

Proposition 3.1. For all m £ {0, . . . , iV} we have u™ 6 RTo . 

Proof. If m £ {0, 1} the result holds by definition. If m £ {2, . . . , TV} we apply the operator div^ to (|3.3p and 
compare with (|3.4[) . We get div ft u™ = 0. Using definition (|2. 12[) we get u™ S RT . ■ 
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Let us show that equation (|3.3p also has a unique solution. Let qt G P\ c n Lq such that A^qh = 0. According 
to proposition 14.41 we have for all qh G Po 

-(A h q h ,q h ) = -(div ft (V hQh),qh) = {^hqh^hQh) = \^hqh\ 2 - 

Therefore we have WhQh = 0, so that qh = since qh 6 ij. We have thus proved the unicity of a solution for 
(|3.3[) . It is also the case for the associated linear system. It implies that this linear system has indeed a solution. 
Hence it is also the case for equation (|3.3[) . Note finally that since u™ G Po n RTo, we have divu™ = for all 
m G {0, . . . , N}. Hence the incompressibility condition (|1.2|) is fulfilled. 

4. Properties of the discrete operators 

We show that the differential operators in (|l.ip - (|1.2p and the operators defined in section [2~3l share similar 
properties. 

4.1. Properties of the discrete convective term 

We define b : H 1 x H 1 — > L 2 . For all u G H 1 and v = (v 1 ,v 2 ) G H 1 we set b(u. v) = (div(«i u), div(w 2 u)) . 
We show that the operator b^ is a consistent approximation of b. 

Proposition 4.1. There exists a constant C > such that for all v G H 2 and all u G H 2 n Hq satisfying 
div u = 

||np b(u,v)-b h (nHToU,np v)||_i, h <cft||u|| 2 ||v||i. 

Proof. We set Uh = Hrt u an d v h = IIp v. Let K G Th- According to the divergence formula and (|2.6p we 
have 



n Po b(u,v)| K = Y Jv(u-n) 



da. 



On the other hand, let us rewrite b^(u^,Vh). Let a G £k H Setting 



VK.L a 



v K si (m • n K ,<r)o- > 

VL„ Si (Uft • n K ,a)tr < 



one checks that v K (u h ■ n. Kf(T )+ + w La (u h ■ n K ^) a = v K ,L a (u/, • tijf^),. Using (|2.1ip . we deduce from (|2. 13|) 
that 



b/,(uh, v;,)|,k = t4t E / V ^ L - ( Uh ' dcr - 
Thus 

n Po b(u,v) -b h (u h ,v h ))\ K = j=r Y / ( v ~ v k,l J (uh ■ n) da. 



Let ip h £P . We have 



(n Po b(u,v) - b h (u h , v h ),tph) = E ^ E (v-VK,L a )(u h -n)da 
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Let a G £™'. We consider the quadrilateral D a denned by x^, XL a and the vertex of a. We set 
Using a Taylor expansion and a density argument (see [14]) one checks that 



D a n K si (u h ■ n K . (y ) a > 
D a n L a si (u h ■ n K a )„ < 



|v - v Ka T a \d<j <Ch 



(L |w(y)|2dy ) 



1/2 



Thus 



(n Po b(u,v) - bft(n RTo u,np v),^) 

<Ch\\u\\ H 2 I \*l>L a -1p 



( E / |Vv(y)| 2 dy] 



1/2 



so that (n Po b(u,v) - b/^IlRToU, n Po v),-0, 1 ) < Ch ||u|| H 2 ||^hl|i,ft ll v l|i- Using then definition ([2T5]), we get 
the result. ■ 

Let v G L°° n H 1 and u G H 1 with div u > a.e. in f2. Integrating by parts one checks that J Q v • b(u, v) efoc = 

i i 2 

In 2~ divuefoc > 0. The operator b^ shares a similar property. 

Proposition 4.2. Let G RTo such that divu/j > 0. For all v/j G Po we Ziawe 

b/,(u ft ,v ft , v ft ) > 0. 

Proof. Remember that for all edges cr G f™*, two triangles K a ct L CT share cr as an edge. We denote by K a 
the one such that u CT • a > 0. Using the algebraic identity 2 a (a — b) = a 2 — b 2 + (a — b) 2 we deduce from 

(EH 



2b h (u h , v h ,v h ) = 2 ^ \a\ w Ka ■ (v K „ - v L J (u h ■ n K<j:<J ) 

= E H(|v^| 2 -|v L J 2 + |v^-v L J 2 )K.n^,.) 



so that 2b fl (u fl ,v ft ,v, 1 ) > J2ae£ l h nt (Wvl 2 - |v L J 2 ) (u h ■ n K ^ a ). This sum can be written as a sum over 
the triangles of the mesh. We get 

2bh(u h ,Vh,Vh) > E l Vi ^| 2 E W\( u h-n K ^,a)- 

KeT h a££ K n£' h nt 

Using finally the divergence formula we get 

2b h (u h ,v h ,v h ) > E 1-^1 \ v k\ 2 / divu^dx>0. ■ 

,7^r Jk 



KeT h 



The following result states that the operator bh is stable for suitable norms. 
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Proposition 4.3. There exists a constant C > such that for all G Po, G Po, u/j G Po satisfying 
diviih = 

\b h (u h ,v h , v h )| < C \u h \ \\vh\\h \Wh\\h- 
PROOF. For all triangle K G T h and all edge a G £k H we have 

(u ft • nif i5 )+ v K + (u/j • njf^); v Lct = (u ft • n K ^) a w K - \ (u h ■ n K ^) a \ (v Lff - w K ). 

Using this splitting, we deduce from (|2.14[) b^u^, v/j, w^) = Si +S2 with 

51 = E Vif-wx ^ |cr| (u/j ■ n^),, 

KeT h <j££ K n£* h nt 

52 = - ^2 w A - • ^ |ct| \(u h -rii<r )0 .) CT | (v^ - v K ). 

By writing the sum over the edges as a sum over the triangles we have 

S 2 = ~ ^ M -T^K,a)A (Vi„ - Vif) • (w Lcr -w K ). 

Using the Cauchy-Schwarz inequality we get 

/ \ 1/2 / x 1/2 

Wsr' J Vest* 

Since u/j G RTo we have [5] the inverse inequality h \\uh\\oo < C \ u h\- Using (|2.1|) and (|2.2p we get 

£ l v ^- v A-J 2 <C £ r CT |v^-v^| 2 <CK|| 2 , 

and in a similar way ^ ff££ »t |w L:T - w A ',J 2 < C || w ?1 1|^ . Thus |S 2 < C \u h \ \\vh\\h \\wh\\h- On the other hand, 
according to the divergence formula 



Si = £ \ K \ ( VK ' / divu ' 1 dx = 



K&T h 

By gathering the estimates for Si and S 2 we get the result. 

4.2. Properties of the discrete divergence 

The operators gradient and divergence are adjoint: if q G H 1 , v G H 1 with v • n|go = 0, we get (v, Vg) 
— (q, divv) by integrating by parts. For and divh we state the following. 
Proposition 4.4. For all v/j G Po and q^ G P" c we have: (vh,Vhqh) = — (qh, div/j v;,,). 

Proof. According to (f2~lg)) 

(v/j, V/j^/j) = £ 1^1 v k • V/j^Ia- = £ v ^ ' ( £ kk/il^) nA, CT 
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By writing this sum as a sum over the edges we get 

(v h) VhQh) = - J! \ a \lh{*-a) (V i<r - VjfJ • n ir „, .+ ^2 W\lh{Xv)v K „ ' ™K„,cr- (4.1) 

On the other hand, using a quadrature formula 

-{q h ,6iY h v ft ) = - ^2 — ^2 (div h v, l )(x CT ). 

By writing this sum as a sum over the edges of the mesh we get 

-(flijdivhVji) = - (^y^ + *( x <t) (div/ l v/ l )(x .) - 9h.(x CT ) (div fe v fc )(x CT ). 

Using definition (|2. 12|1 and comparing with (|4.ip we get the result. ■ 

The divergence operator and the spaces Lq, Hq satisfy the following property, called inf-sup (or Babuska-Brezzi) 
condition (see [9] for example). There exists a constant C > such that 

inf sup _%^)>G (4.2) 

geLl\{0} v£Hj\{0} IM|l|<7l 

We will now show that the operator div^ and the spaces Po n Lq, Po satisfy an analogous property. The proof 
uses the following lemma. 

Lemma 4.1. There exists a constant C > such that 

Vq h EPrnLl sup - i9H f V ^ h) >Ch\\q h \\ lth . 

v h eP \{o} \\Vh\\h 

PROOF. If q h = the result is trivial. Let q h e P™ c n Lo\{0}- Let v h = V h q h E P \{0}. Using proposition WJ\ 
we have 

-(q h ,div h v h ) = (v/,, V h qh) = \V h qh\ 2 = \VhQh\ Wh\- 
Using and we get -(q h , div h v h ) > C h\\q h \\i,h \Wh\\h- ■ 
We now state the result. 

Proposition 4.5. There exists a constant C > such that for all qh £ P\ c fl Lq 

(q h ,div h v h ) 
SU P ii ii >C\Qh\- 

VhGP \{0} \\Vh\\h 

PROOF. If qh = the result is trivial. Let q h e P" c D Lq\{0}. According to (I4.2|) there exists v G Hq such that 

divv = -% and ||v||i<C|g h |. (4.3) 
We set = ilpjv. We want to estimate — (g/ l ,div/ l (IIp v/ l )). Since Vhqh £ Po we deduce from proposition 

-(g/i,divh(IIp v/i)) = (Jl Po v h , Vhqh) = {^h,V h q h ). 
By splitting the last term we get 

- (q h , div h (U Po v h )) = (v, V h q h ) - (v - v h , V h q h ). (4.4) 
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We bound the right-hand side of Using (|2.7p and (|4.3p we have 

|v - v h | = |v - n P jv| < Ch ||v||i < C ft 
Thus, using the Cauchy-Schwarz inequality, we get 

Kv-Vfc.Vfc©,)! < C/i|%| |V fc g fc | < Cfc|©,| ||«fc||i,h. 
We estimate the other term as follows. Integrating by parts we get 

(v ; Vfegfe) = — (g/»divv) + q h (v -n KtCr )da. 

KeT h ae£ K a 

We have — (cfri, divv) = jg^j 2 thanks to (|4.3[) . On the other hand 



X! X! / 9h (v • n K , ff ) da = ^ / g fc (v ■ n^ i<T ) 



KeT h ae£ K a ae£l 

since v|an = 0. Using [2] p. 269 and (|4.3|) we have 



<C'ft||v||i||g fc ||i, h <C , fc|g fc l||«/ fc ||i lh . 



Hence we get (v, V h q h ) > (\qh\ — Ch \qh\i,h) \qh\- Thus we deduce from (|4~4)) 

-(gfc.divfctHpoVfc)) > (|%|-C^||%||i,ft)|%|. (4.5) 

We now introduce the norm \\.\\h- We have v/ t = IIpcv £ PJ C H 1 . From [5] p. 776 we deduce ||IIp Vft||/, < 
C \\vh\\i- Since Upc is stable for the H 1 norm, using (|4.3|) . we get 

|]v ft || a = Hnpjvlli < HvlU < CI^l- 

Therefore ||np Vfc||/, < C | | . Using this inequality in ()4.5|1 we obtain that there exists C\ > and C 2 > 
such that 

-(q h ,diVh(Kp Vh)) > (C\\q h \ - C 2 h\\q h \\ 1 . h ) \\U Po v h \\ h . 

We deduce from this 

(q h ,div h v h ) 

sup n n > Ci\q h \ - C 2 h\\q h \\i, h . 

v h ep \{0} \\vh\\h 

Let us combine this result with lemma I4TT1 Since 



Vt > , max (Ct , d |%| - C 2 t) > |g h | , 

we finally get the result. 

4.3. Properties of the discrete laplacian 

We recall from [14] the coercivity of the laplacian operator. 
Proposition 4.6. For all U/j G Po and Vh G Po we have 

-(A h u h ,u h ) = \\u h \\l, -(A h u h ,v h ) < \\u h \\h\\v h \\h- 
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5. Stability of the scheme 

We first prove an estimate for the computed velocity (theorem 15. ip . We show a similar result for the 
increments in time (lemma l5.2p . Using the inf-sup condition (proposition 14. 5| ) . we infer from it some estimates 
on the pressure (theorem l5.2[) . 

Lemma 5.1. For all m G {0, . . . , N} et n 6 {0, . . . , N} we have 

( u h , ^hPh) = 0, \u h | - |u h I +|u h -u h | =0. 

Proof. First, using propositions 13. II and 14.41 we get (u™, VhP^) = — {p\i div^u™) = 0. Also, we deduce from 

4 A- 

2 (Uh , u h - u h) = — g-( u h i Vft(p h - p h )j = 0. 
Using the algebraic identity 2 a (a — b) = a 2 — b 2 + (a — b) 2 we get 

2 «, u™ - o - |<| 2 - Kl 2 + K - <T = o. 

We introduce the following hypothesis on the initial data. 

(HI) There exists C > such that \u° h \ + |u£| + k\V h pl\ < C. 

Hypothesis (HI) is fulfilled if we set u° = Hrt u o and we use a semi- implicit Euler scheme to compute u^. 
We have the following stability result. 

Theorem 5.1. We assume that the initial values of the scheme fulfill (HI). For all m G {2, . . . , N} we have 

m 

Kf + kY J \\K\\l <c. (5.i) 

n=2 

PROOF. Let m e {2, . . . , N} and n € {1, . . . , m — 1}. Taking the scalar product of (|3.2[) with 4 fc u^ +1 we get 

+4 fc b fe (2 u£ - u', 1 - 1 , u» +1 , < +1 ) + 4 fc (Vfcpjf, u™ +1 ) = 4 fc (f£ +1 , u« +1 ). (5.2) 
First of all, using lemma l5Tl and proceeding as in [10], we get 

4fc (V\ Sflr-^ug + ur 1 ^ = |u „ +1|2 _ |u n |2 + |2u n + l _ u » |2 _ |2< _ 

+ |u» +1 - 2< + u"- 1 ! 2 + 6 |u» +1 - < +1 | 2 . 
According to proposition 14.61 we have — |M| (A^u^ , u£ + ) = ^ ||u^ +1 ||^. Also, according to lemma |5~T1 and 

4fc(V^,< +1 ) = 4fc(V^< +1 -< +1 ) 

= ^ (|VK +1 I 2 - IVKI 2 - iv^r 1 - VKI 2 ). 
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Multiplying equation (|3.4p by 4fc V/j(p^ +1 — and using the Young inequality we get 

According to proposition 14.21 we have 4fcbh(2uJ i 1 — , u? + ) > 0. At last using the Cauchy-Schwarz 

inequality, (|2.3[) and (|3.1[) we have 

4fc(fr i ,< +i ) < 4 fc|fr i i < cfc iifn c(0iT;L2) k +i iu. 

Using the Young inequality we get 

4 fc (f™+ \ < 3 fc ||ar 1 1| 2 , + c k ||f ||2 (o T;L2) . 

Thus we deduce from (|5.2|) 

K +1 I 2 - K| 2 + 12< +1 - <| 2 - |2< - ur 1 ! 2 + K+ 1 - 2uj} + ur 1 ! 2 
+3 k +1 - u ;; +1 | 2 + fc K+ 1 !! 2 , + ^ (|v^ +1 | 2 - |v^| 2 ) < c*. 

Summing from n = 1 to m — 1 we have 

m— 1 m — 1 /i/2 

K| 2 + |2ur - u™- 1 ! 2 + 3 £ - < +1 | 2 + fc 2 K +1 l£ + — ivrfl 2 

n— 1 n— 1 

< C + 4 KI 2 + |2uJt - < | 2 + A: 2 1 V^| 2 . 

Using hypothesis (HI) we get (15.11) . ■ 

We now want to estimate the computed pressure. From now on, we make the following hypothesis on the data 

fGC(0,T;L 2 ), f t G i 2 (0,T;L 2 ), u eH 2 nHj, divu = 0. 

One shows that if the data Uo and f fulfill a compatibility condition [11] there exists a solution (u,p) to the 
equations (|l.ip - (|1.2| such that 

ueC(0,T;H 2 ), u t G C(0,T;L 2 ), Vp G C(0, T; L 2 ). 

We introduce the following hypothesis on the initial values of the scheme: there exists a constant C > such 
that 

(H2) {ul-uol + ^H-uWU + lpl-pihXKCh, \ui-ul\ <Ck. 
One checks easily that this hypothesis implies (HI). We have the following result. 

Lemma 5.2. We assume that the initial values of the scheme fulfill (H2). Then there exists a constant C > 
such that for all m G {1, . . . , N} 

Proof. Using proposition 14. II one proceeds as in [14]. The difference lies in the way we bound the term Vhp\- 
We use the splitting 

v\ = (p\ - Upncpih)) + (Tl PrP (ti) - p(h)) + P (h). 
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Using an inverse inequality [2] we have 



|v fc (pI - Up^pih))] <j\pI- nj»cj>(ti)| < j (\pl-p(h)\ + \ P (h) - nj»cp(ii)|) . 

Using (|2.10p and hypothesis (H2) we get 

|V fc (pi - U PrP ( tl )) | < C ||p(ti)||i < C |bllc(o,riHi). 

According to (|2"TT0|1 we also have \Vh(p(h) ~ n p«°p(*i))| < C < C ||p||c(o,T;m)- Lastly |Vp(ti)| < 

lbllc(o,r ; Hi)- Thus we get \V h pi\ <C, u 

Theorem 5.2. We assume that the initial values of the scheme fulfull (H2). There exists a constant C > 
such that for all m G {2, . . . , N} 



n=2 



< c. 



Proof. Let m € {2, . . . , A}. We set n = m — 1. Using the inf-sup condition (|4. 5[) and proposition 14.41 we get 
that there exists 6 Pq\{0} such that 



C\\v h \\ h \p n h +1 \ < -( P l +1 ,dW h v h ) = (V h p n h + \v h ). 
Plugging p.4p into (|3.2p we have 



(5.3) 



n+l 



3n n+1 - 4n n 4- n" _1 1 ~ ~ 

Re h ~ ^ h ~ h ' h 



n n—1 ~n-\-l\ , rn-\-l 



2k 



so that 



3u 



n+l 



2k 



Re 



(V^ +1 ,v,) = - 

-b,(2<-ur i ,u;: +i ,v, l ) + (f^+ i ,v,) 

Thanks to proposition 14.31 and theorem 15. II we have 



|b,(2<-<- 1 ,u™+ 1 ,v /l )| < (2|<| + |u;r 1 l) IIu^IUKII^ciIu^iuikiu. 

According to proposition 14.61 we have ^A/jU^ +1 , VfA < ||u^ +1 ||h ||vft||ft. Using the Cauchy-Schwarz inequality, 
(j2~3"]) and ([57T]) we have 

and in a similar way 



(^ X ^h)<K +L \Wh\<C\v h \<C\\v h \\ h 



3<+ 1 -4u^ + u»- 
2fc 

Thus we get 

(V hP l+\v h ) <C + C 
By comparing with (|5 .3[) we get 



< C 



3u" +1 ~4u» + ur 1 



2A- 



\\Vh\\h- 



|3< +1 -4ug + u»- 
2fc 



+ K +i IU KIU 
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Squaring and summing from n = 1 to m — 1 we obtain 



<C + Ck 



m—1 



3u! 



-4uJ 



u 



ra-1 |2 



4fc 2 



Ck 



The last term on the right-hand side is bounded, thanks to theorem 15. II And since 

3< +1 -4< + uJT 1 - 3« +1 - <) - « - u"- 1 ) = 3«K +1 - <K 



we deduce from lemma 15.21 

jo n+l a n , „n-l|2 m l^.nia 

' 4fc 2 ~~ ^ fc 2 ~ 

n— 1 n— 1 
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